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Abstract 

In this paper perturbative unitarity constraints for the parameter spaces of two variants of Two 
Higgs Doublet Model with a Z2-symmetric potential are reconsidered, including the latest LHC data 
on existence of a 125 GeV Higgs boson. The analysis does not depend on the Yukawa interactions. 
We focus on the model in which only one of the doublets develops a non-zero vacuum expectation 
value (vev) - the Inert Doublet Model (IDM) and the Mixed Model where both of the doublets have 
non-zero vevs. Positivity constraints, perturbative unitarity condition and the constraints following 
from electroweak precision tests together with conditions determining type of the vacuum are included 
in the analysis. Bounds on quartic parameters of the potential, which are universal for both of the 
models, are given. For the IDM constraints on quartic couplings and masses of scalars as well as their 
implications for Dark Matter scenarios are presented. In the Mixed Model a bound on the value of 
tan/3, 0.18 < tan/3 < 5.59, is found. It depends on the mass of the Higgs boson and is independent 
of the Yukawa sector. In all the cases overall numerical bounds and allowed regions in the parameter 
spaces are presented. 



1 Introduction 

Once the most general ^-symmetric potential of a Two Higgs Doublet Model (2HDM) is written down, 
still different physical models are possible, for a recent review see [I]. Among these are the Inert Doublet 
Model (IDM) [2j|4] in which only one of the doublets develops non-zero vacuum expectation value (vev) 
and the Mixed Model where both of the doublets have non-zero vevs. 

The aim of this work is to constrain the parameter spaces of these two 2HDMs in a consistent way, 
using the following: perturbative unitarity condition, conditions determining the type of the vacuum 
(determining the validity regions of the Mixed Model and of the IDM) , vacuum stability conditions and 
the Electroweak Precision Tests (EWPT) with use of the S and T parameters. Recently a scalar particle 
of mass around 125 GeV has been discovered at the LHC (5j|6J . We assume that the only Higgs boson of 
the IDM corresponds the discovered scalar particle, while in the Mixed Model we consider two candidates 
which may play the role of the Standard Model (SM)-like Higgs boson. The analysis is limited to the 
scalar sector and does not involve constraints from the Yukawa interactions. 

The perturbative unitarity condition was explored in the Mixed Model, following the approach of 



Refs. [7j|9], by many authors, see Refs. 10-20 . Here we improve those analyses by explicitly including the 
conditions determining the existence of the Mixed vacuum and by considering the full tree-level high energy 
scattering matrix. An extensive analysis of the parameter space of the 2HDM with two non- vanishing vevs 
was performed in Ref. |21| . However, this was done under certain assumptions, namely: soft Z2 violation, 
degeneracy of masses of A and H and the SM-like scenario. Here we consider a Z2-symmetric model, thus 
the results of [2l] can not be directly transferred to our caseQ Interesting analysis of the parameter space 



of 2HDM (II) has recently appeared, 22 . The consequences of the stability and perturbativity conditions, 
that are assumed to be valid up to a certain cut-off scale, are analyzed there. Also the experimental data 
is incorporated. However, the main focus of |22| is on the dependence of the allowed parameter space 
on the cut-off scale, which is different than in our case. In Ref. |23| oblique parameters in the most 
general CP-violating 2HDM (Mixed) were studied, their possible values were discussed (also higher order 
parameters V, W and X were analyzed) and some bounds on the mass of the charged Higgs boson were 



1 In Ref. |2l| the Z2 violating parameter is fixed to a nonzero value so it is not possible to consider the Z2 conserving 
limit. 
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found. We focused more on the impact of the EWPT on the allowed regions of all masses in the 2HDM 
and on constraints for tan (3. 



The unitarity constraints for the IDM were studied by us 24 
Refs 



26 27 . The EWPT were analyzed for the IDM in Refs. 



25 and were included also in analyses of 
|3[ |27[|28] . Ref. (27) combines a wide ranee 



of constraints for the IDM, including various experimental results. However, there the main focus is on 
the possibility of accommodating in the IDM a heavy Higgs boson. Here, as was mentioned before, we 
will focus on a 125 GeV Higgs boson. 

The paper is organized as follows: in section [2] the model is briefly introduced, IDM and Mixed 
Model are defined and constraints relevant for further analysis are presented. In section [3] the standard 
perturbative unitarity approach as well as the oblique parameters are revised and a short description of the 
method used to obtain the results is given. Next sections contain the results of the analysis: in section [4] 
bounds on the quartic parameters of the potential are given, in sections [5] and [6] constraints for IDM and 
Mixed Model are presented. Section [7] briefly summarizes the obtained results. 



2 The Models 
2.1 Potential 

We consider a 2HDM with the following potential: 



V 



+A 3 (4&0(4>fo») + H<t>Un)(<t>Us) + |A 5 \(4>Ud) 2 + (4}<f>s) 



(1) 



The parameters m^ 1; m| 2 and Ai . 
that A 5 < (29|{32] 



. A4 are real numbers and without loss of generality we can assume 



Note that this potential is symmetric under two Z 2 discrete symmetries, the so-called D-symmetry 
under action of which: cf>s — > cj>s and <f>D — > —4>d and S'-symmetry which acts as follows: (f>s — > —<l>s and 
4>d — > 4>d [29] . The Z2 symmetry prevents the existence of Flavor Changing Neutral Currents (FCNC) at 
the tree level. 

The requirement that the potential is bounded from below, which is necessary for a stable vacuum 
state to exist, leads to the following conditions (vacuum stability conditions/positivity constraints) [2j: 



Ai > 0, A 2 > 0, A 3 + VA1A2 > 0, A345 + VA1A2 > 0, 



(2) 



where A345 = A3 + A 4 + A 5 . 

A state of the lowest energy (for which the potential approaches its global minimum) is called a vacuum 
state. The potential Q can develop five different types of minima. It depends on the values of parameters, 
which of the minima is the global one. Models in which scalars interact according to potential ([lj but 
distinct vacua are realized, differ significantly in physical content as they develop different particle spectra. 
Here we consider two of themj^j to be presented below. 



2.2 The Inert Doublet Model 

The Inert Doublet Model (IDM) is defined as a 2HDM with the potential (JT|) , an Inert vacuum state and 
Model I of Yukawa interactions |2}j4]. According to the Model I, only <fis couples to fermions. We fix 
model of Yukawa interactions for completeness, our results are independent of Yukawa sector (only DM 
considerations need definite interactions with fermions). 

The Inert vacuum is realized when only one of the doublets has a non-vanishing vev. We take (</>s) = 

9. V (4>d) = 0, consistently with Model I. The doublets can be decomposed as follows: 



<t>s 



<Pl 




2 In the remaining ones unphysical effects, like massive photons and massless fermions, may be present [29]. 



(3) 



In this case the matrix of the second derivatives of the potential is diagonal, so the fields G , H , h, G 
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H, A are mass eigenstatesj^] Their masses read: 

Ml = m\ x = X lV 2 , 
M% ± = i(A 3 « 2 -mi 2 ), 

Ml = 1(^345^ - "4) = + § A 45^, 

M| = i(A 345 v 2 - m| 2 ) = M| ± + iA 45 « 2 , 

where A^" 45 = A3 + A 4 — A5, A 45 = A4 — A5, A45 = A4 + A 5 . G* and G are massless Goldstone bosons. 
The parameters Ai and A3 . . . A5 can be easily expressed in terms of masses of physical particles and the 
parameter m\ 2 - 



The Inert vacuum can be realized only if the following conditions are fulfilled |24 29 



Ml M%, Ml M£ ± > 0, ^>^, (5) 



where the masses of the particles are defined by Eq. Q . 

In the IDM the unique Higgs boson is the h particle. In accord with the Model I of Yukawa interactions 
it couples (at the tree-level) to fermions just like the SM Higgs and so does to the gauge bosons. Thus, h is 
SM-like and we assume that it corresponds to the boson discovered at the LHC 5|6 and set Mh = 125 GeV 
data. 

As the Inert vacuum preserves the D-symmetry and so does the Lagrangian, this symmetry is exact in 
the IDM. Therefore the lightest D-odd particle is stable and can be a good dark matter (DM) candidate 
given that it is electrically neutral. Here we assume that H is the DM candidate [29], so the condition 
A 4 + A 5 < must holdQsee Eq. Q. 

2.3 The Mixed Model 

The Mixed Model is defined as a 2HDM with potential described by Eq. ([T|, vacuum state of the Mixed 
form and Model II of Yukawa interactions. We specify the form of interaction with fermions only for 
completeness, it does not affect our analysis. 

The Mixed vacuum state is realized when both of the doublets develop non-zero vacuum expectation 

values (vev's): (<f> s ) = (<M = ( ^ ), v s , v D ^ 0, v 2 = u| +v% = 246 GeV. The fields can be 

represented as follows: 



p + s ) fo= Pi 

-^{vs + Ps + ixs) J ' \^(vd + Pd + iXD) 

± 



(0) 



Physical particles (mass eigenstates) arc mixtures of the p K , pk and xk ( see f° r example 29 -31 ), namely: 



where R denotes a rotation matrix of an angle a or ft a £ (— 7r/2, 7t/2), f3 € (0, 7r/2) and tan/3 
G and G stand for the Goldstone bosons and H ± , A, H, h are the physical particles of masses: 

M l± = -H^ + x 5 )v 2 



(7) 



M\ 



M h = T T+i^p ( A i + A 2 tan 2 (3 + ^(X 1 - A 2 tan 2 /3) 2 + 4A 2 45 tan 2 /?) , 
M h = ii+i^p( X i + A 2 tan 2 (3 - V(Ai-A 2 tan 2 /3) 2 + 4A| 45 tan 2 fj . 

The relative (with respect to SM) strength of the coupling of h and H to the gauge bosons is controlled 
by (3 — a. Namely, sin(/3 — a) corresponds to the hW + W~ and hZZ vertices and cos(/3 — a) to HW + W~ 
and HZZ vertices. The role of a SM-like Higgs particle can be played both by h (when sin(/3 — a) = 0) 
and by H (when cos(/3 — a) =0). 



3 The particles having their origin in the </)£) doublet {H~, H, A) are sometimes called the dark (inert) scalars. 
4 This condition holds also in the Mixed Model, see Eq. Then, the last two positivity conditions in Eq. jSJ, reduce to 



the last one. 
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The global minimum of the potential is of the Mixed type (i.e. Mixed vacuum is realized) only if the 

A345 m ll 



following conditions are satisfied |24 29 

m\ x \ 2 - A 345 m| 2 



4 



X 2 

A 345 



> o, 



A1A2 

A 4 + A 5 < 0, A 5 < 0, AxA 2 



■ D 



\ 2 

A 345 



AiA 2 - A§ 45 
> 0. 



>0, 



(8) 



3 Method of the analysis 
3.1 Unitarity constraints 

In analyzing the unitarity constraints we follow the standard approach of Refs. [7} |l0{ [l2||20] . From 
unitarity of the S matrix it follows that the full partial wave amplitudes should lie on the so-called 
Argand circle. In particular, in the limit of energy transfer much grater than the masses of the particles 
involved in the scattering, inequality |3ft(a^(s))| ^ \ holds for elastic scattering partial wave amplitudes 
a"'. For the theory to be perturbative, the zeroth order amplitudes should not lie too far from the circle, 
which can be assured by assuming that 



15-17 20 



|R(4 0) (*))| < 



(9) 



where a^\s) denotes the tree-level amplitude of the s-wavej^] 

It is well known that a possible threat to unitarity is due to the scattering of the longitudinally 
polarized vector bosons. The scalar particles are supposed to unitarize the scattering amplitudes. Due to 
the equivalence theorem 1 8 , 9 33 , in the high energy limit it is sufficient to take into account scattering 
of the Goldstone bosons in place of the longitudinally polarized vector bosons. In addition, in the same 
regime, the terms corresponding to cubic couplings are suppressed and only the quartic terms are relevant. 
Moreover, as the physical fields (mass eigenstates) are obtained from the original fields appearing in the 
basic Lagrangian by means of a unitary transformation and only the eigenvalues of the scattering matrix 
are i mpo rtant for the following analysis, it suffices to consider the scattering matrix between the original 
fields 



15 



. So, using the high-energy formula for (s), the unitarity condition, Eq. (JoJ) , can be reexpressed 



in terms of the eigenvalues Aj of the scattering matrix: 

I A, I s$ 8tt. 



(10) 



We consider the full tree-level high energy scattering matrix (of dimension 25) of the scalar sector 
There are 14 neutral channels 



overlooked before 



24 



2.-, 
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charged channels |16||T7] and 3 doubly charged channels, which were 
j Diagonalization of the scattering matrix leads to 12 distinct eigenvalues being 
functions of the parameters A^. As the doubly charged channels do not introduce new eigenvalues, they 
are the same as those from Refs [15-17 20 . Applying the perturbative unitarity condition (10) to these 
eigenvalues yields a set of inequalities for A, . 

Equivalently, different set of parameters can be chosen. The IDM can be parametrized by the parame- 
ters Xi and m| 2 (m^ is just a function of Ai, Eq. Q) or, for example, by the masses of the scalar particles, 
A2 and m22- The Mixed Model can be described by Ai and tan/3 or equivalently, by the masses of the 
scalars, tan/3 and sin a (or sin(/3 — a)). Rewriting the inequalities following from the unitarity condition 
in terms of the latter sets of parameters gives bounds on the masses of the scalars and, possibly, on the 
other parameters. 



3.2 Oblique parameters 

With use of the oblique parameters the contributions from the New Physics (NP) to the electroweak 
processes can be easily tracked. We follow the definitions from [3~H|35] (see references therein). Namely: 



T 



lfA ww (0) A zz (0) 
a 



M w 



As 2 w c 2 w ( A zz (M 2 )-A zz (0) dA^q 2 ) 



a 



Mf 



dq 2 



r 2 



dA lZ {q 2 ) 



cw s w 



dq 2 



q 2 =0 



5 Also a different form of the perturbative unitarity condition was used in the literature, namely |aQ°^| ^ 1 [8 14fll8[|l9 
This of course introduces some differences between results of different papers. 



fn Refs. 18|l9 different classification of states was used, also scatterings between states of hypercharge 2 were considered. 
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where a — e 2 /(47r) is the fine-structure constant, sw = sin 9w, c-w — cos^vk are the sine and cosine of 
weak mixing angle and Ayv is defined as follows: 

n^ v , = g^A vv ,{q 2 ) + q»q»B vv ,{q 2 ). 

Here H^yi is the vacuum-polarization tensor and VV denotes the divectors: 77, jZ, ZZ or WW. 
Moreover, one has to remember that Avv'ist 2 ) contains only contributions from the NP, namely: 

A vv ,{q 2 ) = A™,{q 2 )-A v ™ l [q 2 ), 

where Ay V ,(q 2 ) denotes the quantity calculated in considered model (in this case 2HDM) and Ay^ v ,(q 2 ) 
denotes the same quantity computed in the SM. 

The values of S and T parameters were found to have the following values |36| 7 (for U fixed to be 
equal 0): 



S = 0.03 ±0.09 
T = 0.07 ±0.08 



(11) 



The correlation of the fit is equal to 87%. Using the programe (38] modified by us we generate the ler and 
2a ellipses. We assume that a point in the {S,T) space fulfills the experimental requirements if it falls 
inside the 2a ellipse. 



3.2.1 Oblique parameters in the IDM 

For the Inert Model the expressions for S and T we have taken from [3] (see also 26 35| j 



T 



F(M 2 H± , M 2 A ) + F{M 2 H± , Mjj) - F{M A , M 2 H ) 



8tt 



log 



M, 



h ref 



(12) 



1 (l M|(Mj-3M|) (Mj\ 1 (M* 
2n\6 M 2 A ~M 2 H g {M 2 H J + 6 g \M 2 H± J 36 ' 3 {M\ 



M 2 H M\ 



(13) 



Mh Tci denotes the reference value of the SM Higgs mass used in the fit of the experimental data (Mh ri 
117GeV). Moreover, function F is defined as follows: 



F(x,y) 



2 



log - for x + y, 

x—y y ' 

for x = y, 



3.2.2 Oblique parameters in the Mixed Model 



For the Mixed Model, using formulas from |34p5] adapted for a two-doublet case, the following expressions 
for the S and T parameters are obtained (compare also |39|): 



T 



1 



167r 2 u 2 a I 



{F{M 2 H± , M\) + sin 2 (/3 - a)F(M 2 H± , M 2 ) + cos 2 (/3 - a)F(M 2 H± , M 2 )+ 
- sin 2 (/3 - a)F{M 2 A ,M 2 H ) - cos 2 (/3 - a)F{M 2 A ,Ml)+ 
os 2 (/? - a) (F(M 2 , M 2 H ) - F{M 2 V , M%)) + 

sin 2 (/3 - a)(F(M 2 ,M 2 ) - F(M 2 ,,M 2 )) - (F(M 2 ,M 2 r J - F(M 2 V , M 2 rcf ))J }, 



(14) 



S = ^ {(2s 2 w - 1) 2 G(M 2 H± , M 2 H± , M 2 ) + sin 2 (/3 - a)G(M 2 A , M 2 H , M 2 ) + 

+ cos 2 (p - a)G(M A , M 2 ,M 2 ) - 2 log M 2 H± + log M 2 A + log M 2 + log M 2 - log M 2 r=f + (15) 
+ cos 2 (/? - a)G(M 2 H , M 2 Z ) + sin 2 (/3 - a)G(M 2 , M 2 ) - G(M htet , M 2 )}. 



7 The values of S and T have been updated recently, see |37| , however, the change hardly influences our results. 
8 Note that in the formulas in |26|34| terms containing and do not appear. They are assumed to cancel, however 

it is not the case unless Mh = Mj, " , which is not always true. 
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The following definitions were used: 

16 5(x + y) 2(x-y) 2 3 f x 2 + y 2 x 2 -y 2 (x-y) 3 \ x 

G(x, y,z)= - — + + + log - + 

6 z z z z\ x — y z 6z z J y 

z 2 — 2z(x + y) + (x — y) 2 9 „ , . 

+ V ^— ^ —f(x + y-z,z 2 - 2z(x + y) + (x- y) 2 ), 



for y > 0, 



f(x,y) = \o for y = 0, 

2^J—y arctan for y < 0, 

A/ x 79 z x 2 / „ x x 2 x 2 „x + y\, x ( n x x 2 \ f(x,x 2 — 4xy) 
G(x,y) = -— + 9--2 — + -10 + 18--6 — + — -9 y - log-+ 12-4- ] 



3 y y 2 \ y y 2 y 3 x-yj y V y y 2 ) y 

3.3 Method of the analysis 

We scanned randomly the parameter spaces of IDM and the Mixed Model checking whether the following 
conditions are fulfilled: 



positivity conditions, Eq. |2]); 
perturbative unitarity condition, Eq. (10 I; 



conditions determining the type of the vacuum: Eq. (|5| or Eq. ^ respectively; 
A 5 < and A 4 + A 5 < 0. 



(16) 



To constrain the masses of the scalar particles we took into account also: 

- EWPT (2a). (17) 

In the numerical analysis the values of parameters Ai were chosen randomly from the following ranges: 
Ai, A 2 e (0, 35], A 5 € [-20, 0), A 4 e [-30, -A 5 ), A 3 e (-VX1X2 - A 4 - A 5 , 35]. 

In the IDM the mass of h is set to Mh — 125 GeV. Remaining masses are in the ranges: Ma, M h ± G 
(0, 1010] GeV, Mh S (0, min(M4, Mjj±)], because H is supposed to be the DM candidate and thus has 
to be the lightest of the dark scalars. A2 is in principle allowed to be in the region (0, 35], however it is 
also constrained by Eq. ^ and (|5| . 

In the Mixed Model the masses are allowed to be in the following ranges M H ± , Mh, Ma G (0, 800] GeV, 
M h G (0, Mh] (see Eq. Q), tan/3 e [0, 60Qand sina e [-1, 1]. 

4 Constraints on the quartic parameters of the potential 

We wanted to find constraints on the quartic parameters of the potential, which would be valid both for 



IDM and Mixed Model, so we imposed on Ai the conditions (16 1 (without the conditions determining the 
type of the vacuum). The results read: 

< Ai < 8.38, 
< A 2 s? 8.38, 

-6.05 < A 3 < 16.44, (18) 
-15.98 < A 4 < 5.93, 
-8.34 < A 5 < 0. 



As the conditions (16) include perturbative unitarity condition, the inequalities above show the region of 
perturbativity of the theory. Most of the bounds are significantly more stringent than the frequently used 
bound: |Ai| < 4ir. 

Some of the parameters appear to be interestingly correlated, examples of such dependencies are shown 
in Fig.JTj On the other hand, some are totally uncorrelated, for example a corresponding plot in (Ai, A2) 
space is a rectangle. 



9 There exists a lower bound on the value of tan/3 (tan/3 ^ 0.29) which comes from the assumption of perturbativity of 
the ibH^ coupling and is valid in the Models I-IV of Yukawa interactions |40| . We do not impose it here in order to observe 
the impact of the considered conditions on the small values of tan 0. 
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A 4 



Vu.i 



A, 



Figure 1: Regions of parameters allowed by constraints (16) (without conditions determining type of the 
vacuum) . 



5 Results for the IDM 

5.1 Constraints on the quartic couplings 

In the IDM the quartic coupling constants are simple combinations of the parameters A^ and can be 
constrained with use of the same conditions as for Ai (A34 = A3 + A4): 



-8.10 s$ A345 12.38, 
-7.76 < A3 45 16.45, 
-8.28 < iA 45 sc 0, 
-7.97 < 6.08, 
-7.77 < A 34 sS 13.25. 



(19) 



The combinations listed above correspond to the couplings between scalars, for example A345 represents 
the coupling of hhHH . A345 also controls the annihilation HH — > h, so the DM relic density is sensitive 
to its value. The correlation between A345 and A2 (e.g. H + H~ HH coupling) is presented in Fig. [lj (right). 

It is important to remember about the correlations between parameters as they may introduce ad- 
ditional constraints. For example, it has been shown recently in Ref. |26] that the enhancement in the 
h — > 77 channel in the IDM is possible only for A3 < (A3 is the coupling between h and H + H~). So, as 
follows from Fig. [T] the enhancement in the h — > 77 channel is only possible for —6 < A4 < 5. 



5.2 Constraints on masses 

In the IDM h plays the role of the SM Higgs boson so we fix its mass to M/j = 125 GeV. 
m| 2 enters the expressions for masses (Eq. Q), but does not enter the scattering matrix (Sec 
the bounds on masses depend on its value. For the case with to 22 = the bounds on the masses of scalars 
following from the constraints ( 16 1 read: 



The parameter 
Thus, 



3.11 



M H ^ 602 GeV, 
M H ± ^ 708 GeV, 
M A < 708 GeV. 



(20) 



We have checked in a large range of m\ 2 , for |to 22 | ^ 10 4 GeV 2 , the results hardly change with respect 
to the case with m| 2 = 0. However, when m| 2 is an extremely big and negative, the allowed regions of 
masses are changed, also lower bounds on masses develop. The regions of masses allowed by conditions 
(16 1 and (17 1 for both of the cases are presented in Fig. [2] in the (Mh±, Ma) and (Mjj±,Mh) planes. 
The graph representing correlations in the (Ma,M#) plane is not displayed as it is very similar to the 
one for (M H ± , Mh) (Fig. [2J right). Note that the regions of masses allowed for the cases with m 22 = 
and m 22 = — 10 6 GeV 2 have empty intersection (see Fi g. |2| right panel). 

It can be seen in Fig. |2]that if apart from conditions ( |16[ ) we also impose the EWPT ( 17 ) , the parameter 
space is further constrained. It is most visible in the (Mh±, Ma) plane. Nevertheless, the overall bounds 
on the values of masses p0| remain unchanged. If we examine the possible values of S and T parameters 
superposed on the la and 2a ellipses, Fig.[3j we see that the EWPT constrain the parameter space 
dramatically. 



7 




1 



— 1 — 1 — 1 — 1 — u Mti+ 

800 1000 n 

Figure 2: Regions of masses allowed in IDM by constraints (16 1 (darker regions) and by (16) together 



with (17 1 (lighter colored regions). In both of the plots the figures in the lower left corner correspond to 

10 6 GeV 2 . 



the case with m 2 ^ 



0, while the figures in the upper right corner correspond to 




Figure 3: Values of the S and T parameters calculated for the values of masses fulfilling constraints ( 16 1 
for m| 2 = (left) and m| 2 = — 10 6 GeV 2 (right). The ellipses correspond to la and 2a bounds. 



5.3 Constraints on the Dark Matter 



IDM was shown [28 41 -43 to accommodate a good DM candidate. In this work we assume that H is the 



lightest D-odd particle and thus is the DM candidate. Only three regions of DM masses are considered 
to be consistent with astrophysical data: M DM < 10 GeV, 40 < Af DM < 80 GeV or Af DM > 500 GeV. 
Although the unitarity bounds on Mjj do not effect DM considerations in the IDM in the low mass 
regimes, they do in the high mass region. If H is in the high mass regime it can get only as heavy as 
602 GeV (Eq. (J20) and Fig.^. Thus the high DM mass scenario is highly constrained. Of course this is 
the case only if we assume that |m 2 2 | is a reasonably small parameter. Otherwise, if m 22 "C 10 4 GeV 2 , 
H can be very heavy, without violation of perturbative unitarity. But then, all the dark scalars are very 
heavy as well (Fig.|2|. 



5.4 



Parameter m| 2 



Although the parameter m| 2 cannot be constrained directly by the unitarity condition, one can derive 
an upper bound on its value, taking into account the condition determining the existence of the Inert 
vacuum. From Eq. (5 1 rr ^= > ^P=, substituting mfj = M 2 , Ai = ^ (Eq. (4)) we arrive at 



< 



'\ 2 M h v. 
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Using upper bound on A2 (Eq. ( 18 1) and values of Mh = 125 GeV and v — 246 GeV we obtain the following 

l 22 : 



ml 2 < 9 • 10 4 GeV 2 



6 Mixed Model 



6.1 Constraints on masses of the scalars 

Perturbative unitary condition for the Mixed Model has been concerned in the past [To}|2"o], both analytical 



(21) 



and numerical bounds on scalar masses have been obtained. Our bounds following from conditions (16) 

Af ff s^ 698GeV, 
M H ± ^ 690 GeV, 
M A ^ 711 GeV, 
M h ^ 499 GeV. 

They are in good agreement with the most precise analytical results of Rcf. |20 [^] Larger differences (up 
to about 12%) appear between our results and results of Ref. [l6|[l7] (s ee also comparison in (20|), which 
are probably due to uncertainites of the numerical method used in |16l|17|. 



On the constraints (21 1, the bounds following from the EWPT can be superposed. The resulting 
regions of masses allowed by the conditions (16) and (17 1 are presented in Fig. [4] It is visible that EWPT 
hardly change the upper bounds (21 1. Only the bound on Mh± is shifted down to 600 GeV. 



In the plot showing values of S and T parameters obtained for the masses allowed by the conditions ( 16 ) 
superposed on the la and 2a ellipses in the (S, T) plane, Fig.[5j it can be easily seen that a vast majority 
of the points in the parameter space allowed by theoretical constraints (16 1 is excluded by EWPT. We 



have checked the possible values of S and T parameters using the program 2hdmc [39 and the obtained 
results stay in agreement with these presented in Fig. [5j We have also checked that adding the conditions 
m| 2 = and M sca i ar ^ Mw to our code brings our results to reasonable agreement with those of 
Ref. [23], which were based on such assumptions. 




M A 



Figure 4: Regions of masses allowed in the Mixed Model by the constraints (16 1 (dark regions) and by 
(16 1 together with (17 1 simultaneously (light colored regions). 



6.2 Constraints on tan/3 



allowed by the conditions (16 1, see 



also 



44 



We also investigated correlations between values of masses and tan / 

Masses of A and H ± do not display any correlation with this parameter, as they do not depend 
Eq. Q. On the other hand, Mh and Mh do exhibit interesting depenence on tan/3. It is 



on tan/ 

especially instructive to observe the correlations of Mh and Mh and tan/3 at various values of sin(/3 — a), 
which controls the strength of the couplings between h and the gauge bosons. To this end we ran our 
program changing the value of sin(/3 — a) by a step of 0.2 from -1 to 1. 

The results for Mh are presented in Fig.[6j It shows the maximal values of Mh allowed by conditions 
(16 1 as a function of tan/3 (so the area below the curve is allowed). Note that, the same curve was 
obtained for each value of sin(/3 — a) that was considered. From Fig. [6] it follows that if we consider a 
particular value of the mass of h (or at least set a lower bound on it), then tan/3 is constrained, both 
from above and below. For example, if we look at Mh = 125 GeV, then it can be read from the graph 



10 Small discrepancies of order 1-3% are due to the uncertainity of the numerical method applied in our analysis. 
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Figure 5: Values of the S and T parameters obtained in the Mixed Model for masses of scalars fulfilling 



conditions (16 1. The ellipses correspond to ler and 2er bounds. 



that 0.5 < tan/3 < 6.5. This bound can be improved by setting Mh = 125 GeV explicitly in the code of 
the program and by imposing also the conditions following from EWPT (17 1. Then the following bound 
is obtained: 

0.18 < tan p < 5.59. (22) 



It should be underlined that this bound is obtained solely from the constraints (16), (17) and the assump- 
tion Mh = 125 GeV, without any assumptions on Yukawa couplings Fj 





Figure 6: Correlation between maximal values of the mass of h and tan/3 in the Mixed Model allowed 
by constraints (16 1. The curve does not change as we vary the value of sin(/3 — a) from -1 to 1. The 



horizontal line corresponds to the mass of 125 GeV. 



We can consider h of the Mixed Model to be a candidate for the boson observed at the LHC. Then it 
is required to be SM-like, so sin(/3 — a) has to be equal (or close to) one. Of course the bound (22 1 holds 



for such case, so we conclude that in the SM-like 2HDM (Mixed), with h playing the role of the SM Higgs 
boson with mass of 125 GeV, tan/3 = is strongly constrained. 

The results for Mh are presented in Fig. [7] Clearly, they strongly differ from the results for M^. Now 
we have different curves for different values of sin(/3 — a). Like in the previous case, each curve represents 
the maximal value of Mh allowed by conditions ( |16| as a function of tan p. The curves for | sin(/3 — a) | ~ 1 
are bell-shaped and similar to the M™ ax (tan/?) curve. However, for sin(/3 — a) — the corresponding 
curve is just a straight line. This means that fixing Mh will not introduce any bounds on tan /3 unless we 
fix sin(/3 — a) to a non-zero value. 

Therefore, if we consider a SM-like 2HDM (Mixed) with H playing the role of the SM Higgs boson, 
with cos(/3 — a) close to one (sin(/3 — a) ~ 0), then tan/3 is not constrained by conditions (16), (17 1 and 



11 In particular, this bound is also true in 2HDM with Model I of Yukawa interactions. 
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Figure 7: Correlation between maximal values of the mass of H and tan /3 in the Mixed Model allowed by 
constraints (16 1. Different curves correspond to different values of sin(/3 — a), varying from -1 to (left 
panel, from the bottom to the top) and from 0.2 to 1 (right panel, from the top to the bottom) in steps 
of 0.2. The solid horizontal line corresponds to the mass of 125 GeV. 



assumption on the mass of H. 

The source of the bound (22 1 in the Mh = 125 GeV scenario and of lack of any bound in the Mjj = 
125 GeV scenario can be traced in the expressions for masses of the scalars, Eq. ^ [44 . It can be 
easily checked that Mh — > as tan/3 — > or tan/3 — > oo for any fixed values of Ai, A2, A345. Thus the 
curve in Fi g. [6] was bound to tend to zero for tan/3 — > and tan/3 — > 00, which was the origin of the 
constraint ( |22[ ). On contrary, this can not be stated about Mh and therefore no bound on tan/3 follows 
from the assumption Mh — 125 GeV. 

7 Summary 

We performed an analysis of the parameter spaces of two different kinds of 2HDMs with a Z2-symmetric 
potential. Using perturbative unitarity condition, vacuum stability conditions and conditions determining 
type of the vacuum as well as EWPT we constrained the possible values of parameters. As a result 
we obtained the regions in the space of quartic parameters where the perturbative description is valid 
and consistent with theoretical assumptions. In the IDM we found the regions in the spaces of quartic 
couplings and masses of the scalars which are allowed by the imposed conditions and we showed that the 
mass parameter m\ 2 has to fulfill the following bound m% 2 < 9 • 10 4 GeV 2 . In the Mixed Model with the 
Higgs boson h of mass 125 GeV we found that tan/3 is strongly constrained, namely 0.18 < tan/3 < 5.59. 
We stress the fact that the latter bound does not require any assumptions about the Yukawa sector. 
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